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PROGRESS  ON  A BASE  FLOW  MODEL  FOR  EXTERNAL  BURNING  PROPULSION 

Warren  C.  Strahle,*  Gopal  Mehta,^  and  James  E.  Hubbartt* 
Georgia  Institute  of  Technology.  Atlanta,  Ga. 


Abstract 

Progress  is  described  on  a model  of  the  turbulent  axi sym- 
metric base  flow  behind  bluff-base  bodies.  The  purpose  of  the 
model  is  to  provide  an  analytical  basis  for  incorporation  of 
external  burning  and  base  injection  effects  upon  the  base  flow. 
An  integral  technique  is  used  for  solution  of  the  problem.  At 
the  current  time,  however,  extreme  difficulties  are  being  en- 
countered in  the  choice  of  field  variable  profiles  to  insert 
in  the  integral  technique.  Several  choices  of  profiles  and 
governing  moment  equations  are  yielding  a singularity  before 
the  rear  stagnation  point  is  reached  in  a downstream  integra- 
tion of  the  equation. 


Introduction 

This  is  an  extension  of  the  treatment  of  Ref.  1 to 
treat  the  axi symmetric  case  with  a more  realistic  base  flow 
and  combustion  theory.  This  extension  is  now  warranted  be- 
cause of  the  experimental  proof,  obtained  elsewhere,  that  the 
concept  of  external  burning  behind  bluff  base  bodies  will  in- 
deed work  in  a propulsion  mode. 

It  is  first  necessary  to  construct  an  axi symmetric,  tur- 
bulent base  flow  theory  without  reaction  which  contains  suf- 
ficient detail  to  eventually  incorporate  the  many  phenomena 
which  will  be  of  interest  in  the  external  burning  case.  First, 
the  reversed  flow  region  must  be  treated  with  enough  detail  to 
include  estimates  of  the  reversed  flow  velocity  distribution 
and  to  allow  for  base  bleed.  The  reactive  bleed  may  be  a 


Paper  presents  progress  made  on  ARO  Contract  Q^tJCQAaZ^sC, 
J2jl3£Lentitled  "Analytical  Investigation  of  External  Burning 
behind  Bluff-Based  Axi symnetri c Bodies." 
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promising  adjunct  to  external  burning  to  improve  performance 
and  residence  times  in  the  near  wake  are  influenced  by  the  re- 
versed flow  velocity.  Second,  several  forms  of  adiabatic  as- 
sumptions should  not  be  used  because  reactive  matter  will  be 
entrained  in  the  near  wake,  and  may  be  introduced  by  bleed. 
Consequently,  the  energy  equation  must  be  used  in  some  form. 
Furthermore,  an  entropy  layer  in  the  inviscid  stream  must  be 
allowed  for,  because  several  concepts  of  external  burning 
would  introduce  the  fuel  by  injection  into  a supersonic  stream, 
causing  injection  shocks.  The  approach  outlined  below  is  ul- 
timately capable  of  treating  these  phenomena. 

Analysis 

The  approach  followed  is  similar  to  that  of  Alber^  but 
allows  for  the  foregoing  complexities  and  the  axisynmetric  ge- 
ometry. The  equations  being  used  are  the  usual  approximate 
boundary  layer  equations  for  turbulent  flow  of  the  near  wake 
with  a turbulent  Prandtl  number  of  unity  and  assuming  a per- 


fect gas.  These  are 

9/8X  (pu)  + 1/r  9/9r  (pur)  =0  (1) 

pu  9u/9x  + pv  9u/9r  = - dp/dx  + 1/r  9/9r  (rr-j-)  (2) 

pu  9H/9X  + pv  9H/9r  = - 1/r  9/9r  (rpj  9H/9r)  (3) 

Tj  = Mj  9u/9r  (4) 

H = h + u^/2  _ (5) 

h = c T p = pRT  (6) 

P 


In  Eqs, (1-6)  p is  the  density,  u is  the  axial  velocity  in  the 
x-di recti  on,  v is  the  radial  velocity  in  the  r-di recti  on,  p is 
the  pressure,  tj  is  the  turbulent  shear  stress,  Uj  is  the  tur- 
bulent "viscosity,"  H is  the  stagnation  enthalpy,  h is  the 
static  enthalpy,  Cp  is  the  specific  heat  at  constant  pressure, 
T is  the  temperature,  and  R is  the  gas  constant. 

An  integral  technique,  a variant  of  the  method  of 
weighted  residuals,^  is  used  to  solve  the  system  of  equations. 
If  we  integrate  Eq.  (2)  from  0 to  6,  where  6 is  the  near  wake 
thickness,  and  use  Eq.  (1), the  familiar  momentum  equation  is 
obtained: 

d/dx  (peUgA2)  + PgUgA^  du^/dx  = 0 
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6 6 

A-j  = (1  - pu/pgUg)  rdr,  Ag  = | pu/p^u^  (1  - u/u^)  rdr  (7) 

0 0 

Here  Ai  is  the  displacement  "area"  and  Ao  is  the  momentum  area. 
Multiplying  Eu.  (2)  by  u and  integrating  across  the  layer, 
the  equation  for  the  first  velocity  moment  of  momentum  becomes 
the  following: 

6 

d/dx  (PgUgAg)  + 2PgUg(A^  - A^^)  du^/dx  = 2 (3u/3r)^  rdr 

6 6 ° 

A3  = f pu/PgUg  (1  - u^/Ug)  rdr,  A^^  = (1  - u/u^)  rdr  (8) 

0 0 

Integrating  the  energy  equation,  Eq.  (3),  and  using  Eq.  (5), 

6 

d/dx  (PgUghgA^)  + 1/2  d/dx  (PgUgA3)  + dHg/dx  purdr  = 0 

6 0 

Ah  = pu/PgUg  (h/hg  - 1)  rdr  (9) 

0 

where  Hg  is  the  stagnation  enthalpy  at  the  edge  of  the  near 
wake  and  dHg/dx  = 0 for  the  present  developments.  Integrating 
Eq.  (1)  over  the  layer 

PgVg6  - d/dx  (PgUgA,)  + 6^/2  d/dx  (pgUg)  = 0 (10) 

Since  dPg/dx  = dp/dx  and  -PgUg  dUg/dx  under  the  approximations 
of  this  analysis  and  the  external  flow  is  isentropic,  pg  = 
pe(ug)  and  hg  = hg(ug). 

In  order  to  use  Eqs.  (7-10),  profiles  of  velocity  and  en- 
thalpy must  be  chosen.  The  profiles  are  chosen  to  yield  alge- 
braic simplicity  while  containing  enough  of  the  essential 
physics  of  the  problem.  To  start  the  integration  from  the 
base,  the  profiles  of  Green^  are  used,  following  Alber.  These 
are 

pu/pgUg  = k-2P  0<r<h 

= 1 - P - P cos  7T(r-h)/)l  h £ r £ 6 
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h/hg  = 1 - 2P^  0 < r < h 

= 1 - - Ph  cos  Tr(r-h)/Jl  h r < 6 (11 ) 

Here  iM  is  the  thickness  of  the  shear  layer  and  h(x)  is  the 
height  of  the  roughly  constant  velocity  region  of  reversed 
flow.  6 = h + «,  and  P(x)  and  Ph(x)  are  profile  parameters. 


The  unknowns  in  Eqs.  (7-11)  are  Si,  h,  P,  P^t  Ug,  and  Vg, 
whereas  there  are  only  four  differential  equations  in  Eqs.  (7- 
10).  Consequently,  two  more  pieces  of  information  are  re- 
quired. One  is  obviously  information  about  the  external  flow. 
This  is  treated  by  Webb's  approximate  method  of  characteris- 
tics,5  which  essentially  ignores  left  running  characteristic 
lines  and  allows  a point-by-point  integration  in  x without  car- 
rying a full  characteristics  set  downstream.  This  method  al- 
lows computation  of  the  exterior  flow  angle  and  the  flow  velo- 
city vector,  which  yields  ve  = Ve(ug).  The  remaining  require- 
ment is  an  additional  differential  equation.  Several  options 
are  oren:  1)  various  forms  of  centerline  equations  as  used  by 
Albert  and  Ohrenberger  and  Baum,®  2)  additional  velocity  mo- 
ments of  the  momentum  or  energy  equations,  3)  other  forms  of 
moment  equations  or  equations  integrated  over  only  part  of  the 
domain.  Since  it  is  hoped  to  carry  the  chosen  profiles  to  the 
point  where  h intersects  the  axis  and  there  is  strong  un- 
certainty concerning  ut  behavior  pointwise,  it  is  considered 
dangerous  in  turbulent  flow  to  use  centerline  equations.  An 
additional  velocity  moment  of  the  momentum  equation  is  also 
considered  dangerous  because  it  would  heavily  weight  the  high- 
speed outer  portion  of  the  flow  field.  This  objection  can 
also  be  raised  against  a moment  of  the  energy  equation,  but  a 
stronger  objection  is  the  algebraic  complexity  introduced.  A 
suggestion  made  by  Peters^  has  therefore  been  followed.  The 
momentum  equation  is  integrated  from  0 to  r = h + i/2.  This 
yields 


d/dx  (PgUgA^)  - u'  d/dx  (PgU^A^)  = p^u^  r'^/2  du^/dx 

+ r'p^au/arL , 


pu^/  gUg  rdr. 


'^5  " 


pu/PgUg  rdr 


(12) 


where  primed  quantities  indicate  values  evaluated  at  r = h + 
1/2,  This  equation  is  reasonably  simple  algebraically  and 
there  is  some  confidence  that  ut  t*®  evaluated  properly  in 
the  high  shear  region  of  the  main  shear  layer. 
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The  turbulent  viscosity  is  taken  from  SchetzS  and  modi- 
fied into  a form  used  by  Alber  with  the  necessary  constant 
being  chosen  to  match  the  eddy  viscosity  to  known  results  for 
an  incompressible  free  shear  layer  under  similarity  conditions^ 
The  result  is 

Py  = 0.04649  PgUg/a  (A2  “ ^u^ 

where  a is  the  base  radius.  This  form  should  be  adequate  un- 
til the  rear  stagnation  point  is  reached,  at  which  point  the 
term  involving  A^  would  be  removed  and  pj  “ A2,  the  momentum 
area,  in  accord  with  Alber 's  work. 

The  differential  equations  of  Eqs.(7-10)  and  (12),  to- 
gether with  the  profiles  of  Eq,  (11),  may  be  manipulated  to 
the  following  system  of  first-order  nonlinear  differential 
equations: 


a.  dyj/dx  = b. 

y2  = h y3  = P y^  = P^  yg  = Ug  (14) 

Guessing  a base  pressure  and  value  of  Pu,  these  equations  may 
be  integrated  downstream  from  the  base  if  the  initial  shear 
layer  thickness  H is  known.  The  other  initial  conditions  are 
P = 1/2,  h = a,  and  Ug,  determined  from  the  guessed  base  pres- 
sure. When  h intersects  the  axis,  this  approach  is  terminated 
for  a more  realistic  set  of  profiles  to  carry  the  computation 
through  the  stagnation  point”  and  critical  point^  singularities. 

Results 

For  checkout  purposes  it  was  first  decided  to  assume  a 
constant  pressure  exterior  flow  with  the  additional  restriction 
that  dUg/dx  = 0.  This  allows  the  dropping  of  one  of  the  dif- 
ferential equations  since  du  /dx  is  now  known.  Furthermore, 
the  approximate  method  of  characteristics  may  be  dropped.  Com- 
putations are  shown  for  a case  in  Fig.  1.  The  assumed  value 
of  the  stagnation  enthalpy  in  the  core  flow  at  the  base  was 
equal  to  the  freestream  value,  and  other  values  are  shown  in 
the  figure  with  pt,  the  assumed  value  of  the  base  pressure.  In- 
deed the  solution  proceeds  smoothly  until  the  shear  layer 
lower  edge  intersects  the  axis.  Notable  is  the  extremely  high 
value  of  the  reversed  flow  velocity  attained,  uq^  % - 2360  fps. 

In  fact,  this  is  locally  supersonic  with  respect  to  the  local 
speed  of  sound,  which  gives  an  indication  of  trouble  for  the 
over-all  metiiod.  The  chosen  base  pressure  ratio  is  close  to 
that  experimentally  known  for  this  case.'®  Consequently,  this 
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a base  pressure  and  integratinq  downstream 
of  the  base.  Constant  pressure  case. 

case  was  not  expected  to  yield  a character  far  from  the  actual 
solution,  at  least  for  a short  distance  downstream. 

To  see  the  effects  of  interaction  of  the  viscous  flow  with 
the  external  flow,  the  full  set  of  equations  was  integrated  to- 
gether with  the  approximate  method  of  characteristics.  Under 
several  different  forms  of  the  fifth  equation,  including  cen- 
terline equations,  the  solution  would  start  as  shown  in  Fig.  1, 
but  a saddle  point  singularity  was  encountered  before  the  shear 
layer  intersected  the  axis.  This  also  occurred  under  several 
different  assumptions  for  the  exterior  flow,  including  Prandtl- 
Meyer  flow.  Although  there  appears  no  good  argument  that  a 
critical  point  should  not  appear  until  after  the  rear  stagna- 
tion point,  it  is  believed  that,since  other  solutions'^*"  have 
never  encountered  an  upstream  singularity,  these  simplistic 
profiles  are  at  fault.  Furthermore,  the  appearance  of  a super- 
sonic reversed  flow  in  the  constant  pressure  case  is  not  con- 
sidered tenable.  Although  the  profiles  used  appear  adequate 
for  integration,  say  one  base  radius  downstream,  it  appears 
that, when  one  nears  the  complex  reattachment  region,  they  are 
inadequate.  Consequently,  current  efforts  are  being  directed 
to  conversion  to  a different  profile  type  before  the  shear 
layer  intersects  the  axis. 
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Discussion 

OHRENBERGER  (TRW  Systems,  Redondo  Beach,  California): 

I could  not  agree  more  with  your  comments  in  regard  to 
experience  with  these  difficulties.  What  you  are  going  through 
right  now  is  essentially  what  we  went  through  for  about  a year 
and  a half  of  trying  all  kinds  of  things.  There  are  many  things 
that  we  learned  regarding  very  similar  problems  that  you  are 
encountering.  This  is  not  the  time  or  place  to  go  through 
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those.  One  thing  that  I might  say  in  passing  is  that  we  did 
find  that  the  use  of  the  centerline  momentum  equation  did  not 
work  for  the  very  simple  reason  that  the  integral  method  re- 
quires that  the  pressure  vary  axially  but  not  radially.  In  the 
laminar  case,  especially  near  the  base,  large  radial  pressure 
gradients  truly  do  exist,  and  we  found  that  applying  the  mean 
axial  pressure  gradient  to  the  axis  flow  simply  did  not  work 
at  all.  In  fact  our  solution  actually  behaved  similarly  to 
what  you  experience  with  this  singular  behavior  before  you  get 
to  the  stagnation  point. 

STRAHLE: 

I appreciate  those  comments  and  I understand.  There  is 
an  additional  difficulty  in  the  turbulent  case  about  using  any 
centerline  equation  at  all,  and  that  is  I do  not  know  what  the 
turbulent  exchange  coefficient  is  on  the  centerline.  That  is 
why  I wanted  to  use  an  integral  equation,  but  then  one  does  not 
pick  up  a singularity  and  so  I think  I am  stuck. 

OHRENBERGER: 

It  is  interesting  that  you  have  tried  a moment  equation 
in  lieu  of  centerline  energy  equation  and  you  did  not  experience 
in  that  case  a singularity  akin  to  the  thermal  singularity  that 
we  experienced. 

STRAHLE: 

No,  we  went  right  through  and  you  cannot  see  it  in  the 
equations  either.  The  beauty  of  the  centerline  equations  is 
that  you  can  see  it. 

SEDNEY (Bal listics  Research  Laboratories,  Aberdeen  Proving 
Ground,  Maryland): 

I wanted  to  ask  about  the  characteristics.  You  recall 
Chapman's  work!  many  years  ago.  For  the  axisymmetric  case,  he 
had  to  stop  before  he  got  to  the  axis  and  so  all  his  results  for 
that  calculation  have  a sting. 

STRAHLE: 


I do  not  understand  the  difficulty  because  the  method 
of  characteristics  never  comes  down  to  the  axis.  We  only  use 
it  in  the  inviscid  portion  of  the  flow  and  that  is  off  the  axis. 

^Chapman,  D.,  NACA  Rep.  1051,  1951. 
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SEDNEY: 


Yes,  but  he  found  he  had  to  stop  it  at  a significant 
distance  off  the  axis. otherwise  the  calculation  would  blow  up. 

STRAHLE: 


It  has  never  really  gotten  to  that  small  of  r yet.  The 
viscous  portion  fills  a good  portion  of  the  axis  region  com- 
pared to  the  base  height,  down  to  about  a half  at  a minimum, 
and  we  have  not  experienced  any  trouble  with  the  method  of 
characteristics.  We  define  the  boundary  layer  thickness  of  the 
viscous  zone  thickness, and  we  use  the  boundary  layer  equations 
between  r equals  zero  and  that  particular  thickness.  That  be- 
comes an  unknown  in  the  problem,  but  it  is  only  at  that  exterior 
edge  of  the  viscous  layer  that  we  apply  the  method  of  character- 
istics. 

TANG  (McDonnel 1 Douglas  Astronautics,  Huntington  Beach, 
California): 

I think  you  have  employed  a two-dimensional  velocity 
profile  for  the  axisymmetric  case,  is  that  true? 

STRAHLE: 


That  is  true.  It  is  a simple  profile,  it  has  some 
of  the  characteristics  of  the  flows  that  were  needed, and  I did 
not  see  anything  that  prevented  us  from  using  it  in  an  axisym- 
metric case.  Of  course,  the  values  of  the  profile  parameters 
will  be  different  when  we  march  downstream  in  the  axially  sym- 
metric case  compared  with  the  two-dimensional  case. 

KOOKER (Bal listics  Research  Laboratories,  Aberdeen  Proving 
Ground,  Maryland): 

I realize  that  in  the  example  case  you  have  discussed, 
the  heat  release  is  assumed  to  be  zero.  However,  the  general 
case  as  shown  in  your  initial  slide  might  contain  a substantial 
heat  release  in  the  fuel  entrainment  region.  Since  your 
equations  do  not  account  for  radial  pressure  gradients,  is  it 
possible  that  the  integral  approach  may  encounter  additional 
difficulties  due  to  a lack  of  communication  in  the  radial 
di rection? 

STRAHLE: 

When  we  go  to  the  exterior  part  and  actually  include 
external  burning,  I am  planning  to  continue  to  use  an  integral 
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technique  but  including  the  effect  of  the  radial  pressure 
gradient.  Now  that  may  be  too  ambitious  a statement,  but  I 
am  not  there  yet.  I am  still  trying  to  get  to  the  rear  stagna- 
tion point. 
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metric  base  flow  behind  bluff-base  bodies.  The  purpose  of  the 
model  is  to  provide  an  analytical  basis  for  Incorporation  of 
external  burning  and  base  injection  effects  upon  the  base  flow. 
An  integral  technique  is  used  for  solution  of  the  problem.  At 
the  current  time,  however,  extreme  difficulties  are  being  en- 
countered in  the  choice  of  field  variable  profiles  to  insert 
in  the  integral  technique.  Several  choices  of  profiles  and 
governing  moment  equations  are  yielding  a singularity  before 
the  rear  stagnation  point  is  reached  in  a downstream  integra- 
tion of  the  equation. 
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